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Abstract 

An elastic layer slides on a rigid flat governed by Coulomb's friction 
law. We demonstrate that if the coefficient of friction is high enough, the 
sliding localizes within stick-slip pulses, which transform into opening 
waves propagating at intersonic speed in the direction of sliding or, 
for high Poisson's ratios, at supersonic speed in the opposite direction. 

This sliding mode, characterized by marginal frictional dissipation, and 
similar to carpet fold propagation, may result in inversion of the frictional 
force direction; at longer time intervals the system demonstrates stick-slip 
behavior. The mechanism is described in detail and a parametric study is 
presented. 
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1 Introduction 


Frictional slip on interfaces, which separate solids with different elastic 
properties (referred as bi-material interfaces) controls propagation of inter¬ 
facial mixed mode cracks in composite materials [Coker et al., 2003], energy 
dissipation in brake systems [Anderson and Knapp, 1990, Moirot et al., 2003] 
and fault slip [Scholz, 2002, Rice, 1983, Kawamura et al., 2012]. Regardless 
the huge difference in scales between these examples, frictional sliding 
demonstrates scale invariant properties, which enable to mimic Earth- 
scale phenomena, like earthquakes, in the laboratory [Xia et al., 2004, 
Ben-David et al., 2010]. However, the scaling laws and the coupling 
of different mechanisms involved in frictional slip are not yet fully 
understood and present an active research area in mechanical, physical 
and geophysical communities [Miiser et al., 2001, Carpinteri and Paggi, 2005, 
Urbakh et al., 2004, Kawamura et al., 2012]. Difficulties in studying the 
elastodynamic frictional sliding on bi-material interfaces originates in 
unknown exact interface laws, complex interaction of the surface and the body 
waves as well as non-trivial interplay of the involved time and spatial scales. 
The fact that the slip occurs on the bi-material interface implies that the frictional 
slip locally changes the contact pressure, which is not the case for similar 
materials. For the Coulomb's friction law, this strong coupling between normal 
and tangential tractions results in a mathematically ill-posed problem or 
flutter instability [Renardy, 1992, Martins et al., 1995, Ranjith and Rice, 2001] 
for a wide range of combinations of materials and friction coefficients. For 
example, for material pairs, for which a generalized Rayleigh interface wave 
does not exist, the frictional slip is well-posed only if the coefficient of friction 
is below a certain limit: in the case of contact between an elastic half-space 
and a rigid flat, this limit is a unity [Renardy, 1992]. Above this limit the 
well-posedness depends on the Poisson's ratio v and the coefficient of friction 
/ [Martins et al., 1995]. 


/ / / / H= const | /JJJJJJH 


V Q = const <scq 


Elastic layer ( E,v,p ) 


H 


Frictional interface 


/////// / Rigid flat / // // // // 


Figure 1: Problem set-up: infinite elastic layer of thickness H with imposed 
period A slides on a frictional interface with a rigid flat being sheared at 
constant rate Vq. 
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In this Letter we demonstrate that the mathematical ill-posedness loses 
its destructive effect for a finite system, whose frictional dynamics appears 
to be governed by the structural dynamics. We report accurate numerical 
results showing that the frictional sliding in the "ill-posed" regime propagates 
as a series of stick-slip pulses. These pulses, after had run a certain distance, 
transform into stick-slip-opening waves; both modes of slip propagation imply 
very low frictional dissipation. The mechanism is described in detail and a 
parametric study is presented. 


2 Methods 

We consider (Fig. 1) an infinitely long elastic layer of height H (Young's 
modulus E, Poisson's ratio v, mass density p ) pressed against a rigid flat 
by a constant displacement Uy < 0 applied on the top and resulting in uniform 
interface pressure p. The contact interface is governed by Coulomb's friction 
law with constant coefficient of friction / independent of velocity and slip 
distance. From the physical point of view, this implies that at small scale 
the interface and adjacent materials are assumed to be time independent and 
not affected by generated frictional heat. The solid is sheared at a constant 
velocity Vg <s c s applied on the top surface, where c s is the shear wave celerity. 
We use the finite element method with implicit time integration to solve the 
elastodynamic equation with contact and frictional constraints. 

We solve this problem for an elastic layer of length A. To imitate an 
infinite elastic layer, periodic boundary conditions were used on lateral sides 
(displacement degrees of freedom of corresponding nodes are coupled). The 
simulations were carried out in in-house displacement-based finite element 
software Z-set [Zset, Besson and Foerch, 1997] 1 . To avoid spurious internal 
wave reflections, we used a regular structured finite element mesh of square¬ 
shaped linear elements with reduced integration (one Gauss point per 
element). To verify the mesh convergence, we used finite elements with the 
side/; = H/32,H/64,H/128. ImplicitHilber-Hughes-Taylor [Hilber et al., 1977] 
time integration scheme was used with a moderate dissipation of high 
frequency modes a = 0.1. The method establishes the following general 
time-discrete equations of motion: 

Mu t+At + (1 - a)Ku t+At + aKu 1 = (1 - a)F^ At + aF‘ fM , 

where M is the mass matrix, K is the stiffness matrix, F ext is the vector of 
external forces, u, ii, ii are vectors of displacement, velocity and acceleration, 
respectively; the upper index refers to time, and At is the time step. Plane 
strain formulation with infinitesimal deformations is used. 

We solved elastodynamic equation of motion 

V-£7 = pi i, 

^or a single case, the simulations were run also in commercial finite element software 
ABAQUS, which properly reproduced the results obtained with our software. 
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where a is the Cauchy stress tensor and ii is the acceleration vector. The 
stress-strain relationship is given by Hooke's law o = ( 1+v ^\'_i v ) tr(f)i + 
where tr(£) is the trace of the deformation tensor, and I is the second rank 
identity tensor. Since we consider small strain elasticity, the strain tensor 
is given by £ = l/2(Vu + uV). The boundary conditions are the following: 
on the top of the rectangle y = H, vertical displacement u y = -0.001 H is 
prescribed, which results in contact pressure p = U,,E(1 - v)/[H(l + v)(l - 
2v)], on lateral sides periodic boundary conditions are imposed for every 
couple of corresponding nodes u x (x = 0, y) = u x (x = A, y) and u y (x = 0, y) = 
u y {x = A, y), these conditions are kept constant during the simulation. On 
the contact interface the classic Hertz-Signorini-Moreau contact conditions are 
used [Wriggers, 2006, Yastrebov, 2013]: 

P < 0, Uy > 0, pUy = 0 


and frictional conditions for Coulomb 

friction [Wriggers, 2006, Yastrebov, 2013] 

\o xy \ ~ f\p\ < 0, \u x \ > 0, {\a xy \ - f\p\)u x = 0, 

where a xy is the tangential traction on the contact surface defined as a xy = t e n, 
where n is the unit outward normal and t is the unit tangent vector. The 
frictional contact constraints are taken into account via a direct method, 
which uses local stiffness matrices to resolve the contact constraints within 
a nested loop separated from the global convergence loop, in which contact 
and frictional forces enter as external forces, the method was developed 
in [Francavilla and Zienkiewicz, 1975, Jean, 1995]. 

The elastic layer was first loaded by vertical displacement u y and next the 
top surface was sheared by u x until \T/N\ = 0.95/, where 


A A 



0 0 


are the tangent and the normal forces, respectively; both initial steps are solved 
quasi-statically (without dynamic effects). Next, a dynamic simulation started 
(initial velocity field is zero), the top surface is sheared at a constant rate 
u x = Vo, this rate is chosen to be significantly smaller than the shear wave 
celerity; for most simulations presented in the Letter we used Vo = 10 _5 c s . 

In almost all simulations we used Young's modulus E = 1000 Pa, mass 
density p = 1000 kg/m 3 , H = 25 m, Poisson's ratio was varied in the range 
v £ [0,0.49] and the coefficient of friction was varied in the range / £ [0.5,2.]. 
The time step was chosen as At ~ 0.5m /c s To verify the validity of results, 
various mesh periods were tested A = 2 H, 4 H, 6H, mass density was varied as 
p = 0.1,10,1000 kg/m 3 , the time step was adjusted proportionally to p 1/2 . In all 
cases the slip dynamics retains all its properties. Different shear rates were also 
tested V 0 /c s = 10 _6 ,5-10 -6 ,10 _5 ,5-10“ 5 ,10 -4 , they can change the slip dynamics 
at longer time scales, but for moderate shear rates the initial frictional drop and 
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the formation of slip pulses and opening waves were observed for the most of 
studied parameters. Higher shear rate Vq/c s = 10 -4 for some parameters may 
suppress the inversion of frictional force. Note that the influence of the shear 
rate was not the focus of the present study and requires a detailed analysis. 

2.1 Remark on the choice of the time step and the mesh density 

The time step in implicit time integration can be arbitrary large as the method 
is unconditionally stable, however, too large time steps smooth out rapid wave 
dynamics leading to slow structural dynamic simulation. But if the time step 
is chosen sufficiently small, of the order of stable explicit time step (Courant- 
Friedrichs-Lewy condition), the wave dynamics is finely resolved and for time 
steps A t < lm/c s , the choice of the time step does not affect the solution for 
quasi-stationary regime, where only major slip pulses survive and secondary 
ones die. 

The mesh size is a more critical parameter, but it affects only the short 
transitional period from the start of sliding up to the localization of slip pulses. 
According to [Martins et al., 1995], for high enough friction coefficients, the 
homogeneous slip is unstable and the amplitude of interface displacements 
grows with an exponent proportional to the wavenumber or, in our case, to 
the inverse of mesh size. This destabilization may result in localization of 
stick-slip pulses spaced by a single finite element. However, this instability 
is only observed within the stabilization cycle before the finite size of the 
system manifests itself. When the system starts to act as a whole, only 
major slip pulses localize (their number is determined by the ratio of height 
to period and the slip front celerity), at this stage the mesh size does not 
change solution characteristics, if a sufficient number of elements is used 
enabling to capture all relevant fields' variations. We tested meshes with 
32, 64 and 128 elements in height, and for the stabilized quasi-periodic 
solution we could capture quantitatively equivalent results. To summarize, 
the dependence on mesh and time step exists, it is slightly more subtle than 
in classic dynamic simulations, but this dependence is eliminated in the stable 
regime when the finite size of the system manifests itself. Note, that to get rid 
of any manifestation of the ill-posedness, one may use Prakash-Clifton regu¬ 
larization of frictional law [Prakash and Clifton, 1993, Ranjith and Rice, 2001, 
Cochard and Rice, 2000, Kammer et al., 2014]. However, we believe that 
introducing the regularized friction law at the interface will not affect the 
qualitative results within the stable regime. 


3 Mechanism of frictional sliding 

At reaching the critical shear stress t = f\p\, the frictional interface, initially 
in stick (pinned) state, starts to slip uniformly. But since this uniform slip 
is unstable for high friction coefficients [Martins et al., 1995], after a short 
transitional period, slip localizes within several slip pulses propagating along 
the interface (see an example in Fig. 2 obtained for v = 0.2, / = 1.5, £ = 1000 
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time (s) 


Figure 2: Spatio-temporal maps showing formation and propagation of (a) 
slip pulse (ii x > 0) and (b) opening wave (u y > 0), which result (c) in a 
drop of frictional reaction (ratio of tangential T to normal force N on the 
interface) down to negative values; longitudinal, shear and slip wave celerities 
are denoted ci,c s ,c v , respectively; maps (a,b) and the frictional force (c, red 
line) are shown for Vq/c s = 10 -5 , frictional forces plotted in gray in (c) are for 
V 0 /c s e [10 -6 ; 10 -4 ]. In (d), the ratio T/N for the same simulation is plotted for 
a longer time interval, f = 0 corresponds to T/N ~ 0.95/. 
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Pa, p = 1000 kg/m 3 , H = 25 m, A = 50 m, Vq/c s = 10 -5 ). Material velocity, 
slip velocity and opening profiles at distinct time moments, which correspond 
to times marked with white dashed lines in Fig. 2, are depicted in Fig. 3 
(an animation of the slip onset and formation of slip and opening pulses is 
provided in supplemental material [Sup]). Evolution of frictional force for a 
longer time interval (Fig. 2,d) is similar to the classical stick-slip behavior. Note, 
however, that here this global stick-slip dynamics emerges from propagation 
of rapid slip pulses within time intervals corresponding to frictional drops, 
in other time intervals the entire interface remains in stick state accumulating 
elastic energy. Note also that the frictional drops, except the very first one, 
occur when T/N ~ 1 regardless the fact that the local coefficient of friction 
f = 1.5. 

The observed local propagation mode is referred as train of slip 
pulses [Adams, 2000, Coker et al., 2003, Bui and Oueslati, 2010]; solitary 
pulses in finite systems were also observed [Baumberger et al., 2002]. For 
moderate Poisson's ratios, slip pulses propagate at intersonic celerity c p 
(c s < c p < Ci), where C; is the longitudinal wave celerity. Life period of many 
of these initially formed pulses is short: being filtered by wave dynamics 
and the resonance of the finite size system, their propagation speed and slip 
intensity decrease before they arrest in a manner rather similar to what was 
observed in [Kammer et al., 2014]. However, the major pulses keep going 
through the interface at the constant speed. After a certain propagation 
distance an opening occurs within the slip pulse. This opening implies that 
the interface can slide almost without local slipping and thus without any 
energy dissipation. One may think of a carpet fold analogy: a carpet can be 
moved along a surface by simply propagating a fold, which requires much 
less energy than sliding the entire carpet [Comninou and Dundurs, 1978]. 
However, here the opening starts at location, where the slip velocity reaches 
its maximum (Fig. 3,c-e), thus a marginal frictional dissipation still takes place. 
The opening pulses demonstrate a self-sustained behavior and have a certain 
"inertia", so that their propagation often results in bringing the contact surface 
ahead of the top one. This unconventional mechanism results in inversion 
of the frictional force direction and thus in instantaneous negative apparent 
friction (Fig. 2,c). It is partly similar to the snap-through instability with 
accompanied negative dynamic stiffness [Lakes et al., 2001]. The frictional 
dissipation is marginal in this slip mode, thus the total energy is approximately 
conserved. The elastic energy liberated within the first frictional drop via stress 
waves, is partly restored in the solid being sheared in the opposite direction 
(Fig. 3,f). However, in the particular case of the softly applied shear load via 
a "pusher" [Ben-David et al., 2010, Radiguet et al., 2013], this stored energy 
will be also liberated. Such a huge drop in frictional force down to negative 
values was observed in molecular dynamic simulations [Sorensen et al., 1996, 
Urbakh et al., 2004]; also, within a simple stick-slip model (a rigid block 
pulling by a spring over a frictional surface under velocity weakening friction 
with the static and kinetic coefficients of friction f s > /; c ) the spring force 
may change the sign if /jt < f s / 2. Evidently, for systems governed by 
Coulomb fricton with the coefficient of friction / independent of velocity 
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and slip distance, frictional drop is impossible in the quasi-static case when 
inertial effects are neglected. For elastodynamic case, it was argued that the 
apparent kinetic friction (the ratio of the tangential to the normal force) at 
bi-material interfaces should be smaller than the local one [Adams, 2000] 
or even vanish completely [Adams, 1998, Ranjith and Rice, 2001]. This is 
possible due to the drop of contact pressure within the slip pulse theorized 
by Weertman [Weertman, 1980]. However, for the best of our knowledge, the 
changing of the frictional force direction was not observed at macroscopic 
scale. 


4 Bibliographical remarks 

The opening waves, also known as Schallamach waves, is a familiar phe¬ 
nomenon in rubber friction [Schallamach, 1971, Gent, 1974, Barquins, 1985]: 
when a rigid slab slides over a rubber substrate at moderately high 
velocity, the relative motion between two solids occurs in the form of 
detachment waves, the energy is dissipated mainly in the bulk of the 
rubber. This phenomenon is qualitatively similar to our opening waves, 
however, the physics of the Schallamach waves, which propagate at 
smaller speeds than what is observed here, includes visco-elastic material 
behavior and adhesion on the contact surface, both are not present in 
our model. The possibility of detachment waves was also speculated 
in [Comninou and Dundurs, 1977, Comninou and Dundurs, 1978] either for 
similar materials or for frictionless contacts, however, some inconsistency 
in energy balance in these results was pointed out in [Freund, 1978]. 
Possibility of detachment was discussed in numerous works [Weertman, 1980, 
Adams, 1995], but obtaining rigorous theoretical results is associated with 
great mathematical difficulties. Nevertheless, the detachment waves on 
frictional interfaces, different from Schallamach waves, were induced in 
atomistic simulations [Gerde and Marder, 2001] and, for the first time, were 
obtained in finite element simulations for a set-up with rotational symmetry 
in [Moirot et al., 2003]. 


5 Stick-slip and stick-slip-opening pulses 

Here, we consider in detail the formation of an intersonic slip pulse and the 
mechanism of its self-sustained propagation [Weertman, 1980, Heaton, 1990, 
Adams, 1995]. The strong coupling between normal and shear tractions 
on bi-material interface creates favorable conditions to reach the frictional 
limit within a local zone and make slip occur and propagate in the 
direction of sliding of the softer material [Andrews and Ben-Zion, 1997, 
Cochard and Rice, 2000, Ben-Zion and Huang, 2002]. Since the study of a 
solitary pulse in a periodic system is problematic and the number of stable 
pulses depends on the ratio A/H and on the pulse speed c p , we consider a 
train of stick-slip pulses propagating along the interface at constant speed 
Cp. Let I denote the spacing between pulses and s the width of the pulse. 
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Figure 3: Instantaneous maps of material velocity in the elastic layer, slip 
velocity and opening on the interface are presented for / = 1.5, v = 0.2, 
V 0 /c s = 10“ 5 at (a) t * 787 s, (b) t ~ 840 s, (c) t « 890 s, (d) t ~ 960 s, (e) t ~ 1010 
s, (f) t ^ 1077 s; deformed state is scaled by a factor 150; dashed line shows the 
configuration just before the slip starts. 


Such type of motion was conjectured theoretically [Adams, 1998, Adams, 2000, 
Bui and Oueslati, 2010], observed in experiments [Coker et al., 2003] and 
in simulations [Cochard and Rice, 2000, Coker et al., 2003]. However, in 
theoretical findings too restrictive conditions were imposed: the mean slip 
velocity on the interface was assumed to be equal to the shear velocity Vq. 
This condition is natural if the motion is assumed stationary, it is equivalent 
to the following relation u x = VqI/s, where u x is the mean slip velocity in the 
pulse. Such a stationary solution is observed for low coefficients of friction. 
However, such a stationary train of pulses cannot result in a huge frictional 
drop bringing the system in a global pinned state. Moreover, we observe that 
the contact surface may be brought ahead of the top one, thus the link between 
the local and the global speeds should be stated as an inequality 

ti x > V 0 l/s. 
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Hence, to characterize the train of stick-slip waves, we have four unknowns 
Cp, s, /, u x and one inequality. 

5.1 Pulse celerity 

As observed in simulations, in most cases the slip front propagates in the 
sliding direction (direction of the imposed shear) at intersonic velocity c p (see 
Fig. 2 for an example, and Fig. 4,b and Table 1 for the summary). Thanks 
to the similarity between propagation of the mode-II crack and the slip 
front [Palmer and Rice, 1973, Andrews, 1976, Svetlizky and Fineberg, 2014, 
Kammer et al., 2015], to find the slip pulse celerity, one may resort to the 
help of results obtained for intersonic cracks at the interface between an elastic 
solid and a rigid substrate [Liu et al., 1995]. It was shown that the sign of the 
product between vertical displacement at a certain distance behind the crack 
tip and the contact pressure ahead of the crack tip is determined by the sign of 
(Cp/c s - V2). Since the pressure ahead of the slip pulse is always compressive 
(negative), the vertical displacement is positive if c p < V2c s , thus propagation 
at this speed is favorable for opening. In almost all simulations, in which slip 
and, in subsequence, opening fronts propagate in the direction of sliding, their 
celerities were found to be confined in the interval c s < c p < V2c s . The slip 
front celerity weakly depends on the coefficient of friction, and, in general, 
the most energetic pulses are obtained for friction coefficients, which result in 
c p ~ a/2c s . 

In Table 1 slip-pulse celerities c p are provided for several combinations of 
parameters, which result in opening waves. Celerity of slip fronts is positive 
if they propagate in the direction of sliding, and negative if they propagate 
in the opposite one. The data are obtained for E = 1000 Pa, p = 1000 kg/m 3 
and A = 2H; the average celerity is measured over a time period of 200 s 
using a threshold in the slip velocity ii x , and was also verified geometrically. 
The slip celerity remains almost constant when a major slip pulse is localized 
and does not change when the opening occurs, at least within the precision 
of simulations and measurements. For v = 0 and / = 2, in contrast to all 
other studied combinations of parameters, the slip front propagates in the 
sliding direction at supersonic speed, however, this result does not contradict 
the theoretical basis [Adams, 2000]. More generally, it is usual in non-linear 
dynamic systems that for certain combinations of parameters several stable 
(unstable) solutions coexist, whose realization depends on initial conditions. 
Probably, here we face a similar situation. 

5.2 Pulse period 

A relation can be established between c p and the period length l. Since the slip 
pulse propagates at intersonic celerity, it radiates a shear Mach wave which 
propagates in the bulk inclined at angle 6 S = arccos (c s /c p ) with respect to the 
rigid flat. Being reflected from the top surface the Mach wave returns to the 
frictional interface and might create favorable stress state to initiate or sustain 
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/ 

0 . 

0.1 

Poisson's ratio 
0.2 0.3 

v 

0.4 

0.45 

0.49 

0.50 







-5.250 

0.80 







-5.253 

0.90 






-2.506 


0.95 






-2.523 


1.00 






-2.523 


1.05 




0.851 

0.825 

-2.511 

-5.248 

1.10 

0.861 

0.829 


0.857 




1.15 

0.843 

0.824 

0.880 

0.858 




1.20 

0.861 

0.921 

0.875 



-2.504 

-5.085 

1.30 






0.708 


1.40 






0.708 


1.50 

0.982 


0.882 

0.774 

0.679 

0.656 

-5.399 

1.80 


0.945 

0.831 

0.764 

0.717 



2.00 

1.213 



0.723 

0.666 

0.659 

0.607 

C s 

0.707 

0.674 

0.645 

0.620 

0.598 

0.587 

0.579 

V2c s 

1.000 

0.953 

0.913 

0.877 

0.845 

0.830 

0.819 

Cl 

1.000 

1.011 

1.054 

1.160 

1.464 

1.948 

4.137 


Table 1: Slip front celerity in (m/s) for various combinations of Poisson's ratios 
and coefficients of Coulomb's friction. In the three last lines, c s , V2c s , c; are 
provided. 


another pulse. Thus l = 2Htan(0 s ) = 2H ^jc^/Cg - 1, and if c p ~ V2 c s then 

l ~ 2H. This hypothesis was numerically verified for / = 1.5, v = 0.2 by 
varying the ratio A/H = 2,4,6 resulting in 1,2,3 major slip-opening pulses. 
The last unknown, slip length s, remains undetermined. 

5.3 Summary 

Opening waves were observed for certain combinations of friction coefficients 
and Poisson's ratios (Fig. 4). Such waves are accompanied with frictional drop 
(but it does not necessarily leads to inversion of the frictional force). This result 
holds for all shear velocities in the considered interval Vq/c s e [10 -6 ; 10 -4 ] (see 
Fig. 2,c). For high Poisson's ratios, within a certain interval of coefficients of 
friction, slip and opening pulses propagate at supersonic celerity c p > c; in the 
direction opposite to the sliding; their supersonic celerity is compatible with 
theoretical results [Adams, 2000]. For higher friction coefficients, slip pulses 
propagate again in the direction of sliding at intersonic celerity. Seemingly 
unrealistic, the supersonic slip propagation does not contradict the causality 
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in the system with maximal signal speed equal to the longitudinal wave speed 
C;. It means simply that the slip advancement is not a result of the local 
deformation due to the slip, but due to elastic waves traveling in the elastic 
layer, some experimental and numerical results on supersonic slip can be 
found in [Coker et al., 2003, Coker et al., 2005, Kammer et al., 2012]. The time 
interval of the frictional drop and its value is not affected by the applied shear 
velocity Vq. The time step and mesh density affect the slip dynamics only 
within the transitional phase as discussed in Section 2.1, but do not alter the 
formation of major slip pulse and its transformation into opening wave if the 
time step is small enough to resolve accurately the wave dynamics. 


6 Conclusion 

Experimental verification of the observed phenomena may be carried out 
on a set-up with rotational symmetry [Moirot et al., 2003], which would 
ensure an initial uniform stress state; non-uniform stress distribution, 
usual for finite size contacts, results in triggering frictional slip before the 
local frictional limit / is reached on the global scale. Thus, in finite 
systems the apparent or global coefficient of static friction f* should be 
smaller than the local friction /, which is compatible with experimental 
observations [Ben-David et al., 2010, Ben-David and Fineberg, 2011]. In this 
light, the behavior of frictional system with a uniform stress distribution 
resembles to the process of supercooling: the slip may not start before the 
entire interface reaches the frictional limit, but when it does, the process is 
rapid and energetic. On contrary, in systems with stress heterogeneities, easily 
triggered slip front may propagate along the entire interface resulting in global 
sliding [Ben-David et al., 2010, Kammer et al., 2015]. 

The sliding mechanism discussed in this Letter adds a new piece to the 
puzzle of elastodynamic sliding under Coulomb's friction law in finite-size 
systems, but the full picture is still missing especially in the region of high 
coefficients of friction and Poisson's ratios. The prospective study is planned 
to focus on the dynamics at longer time intervals. 
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Figure 4: (a) Intersonic opening waves propagating in the direction of sliding, 
which are accompanied with high frictional drop, were observed for Poisson's 
ratios and friction coefficients marked with a square (double square for cases 
with inversion of frictional force), diamonds mark supersonic opening fronts 
propagating in the opposite direction, open circles correspond to parameters, 
for which opening was not observed; parameters corresponding to theoretical 
flutter instability for slip without opening [Martins et al., 1995] is located in the 
wedge-zone bounded by a dash-dotted line, (b) shows celerities of opening 
waves for different combinations of Poisson's ratios and friction coefficients; c/, 
c s as well as the critical value [Liu et al., 1995] V2c s are also plotted. Examples 
of slip localization are shown in spatio-temporal maps (black zones correspond 
to stick state, colored zones correspond to slipping) (c-f) for (v = 0, / = 1.5), 
(v = 0.2, / = 1.1), (v = 0.3, / = 1.1) and (v = 0.45, / = 1), respectively; the 
latter combination of parameters results in supersonic slip and opening pulses 
propagating in the direction opposite to sliding. 
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